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ABSTRACT 


High  speed  photographic  coserage  of  the  head  and  neck  articulation  of  a  living  human 
volunteer  both  for  voluntary  motion  and  for  oblique  and  lateral  impact  was  reduced  and 
analysed. 

This  analysis  identifies  a  range  of  four-pivot  mechanical  linkages  that  can  simulate 
the  head  and  neck  kinematics  of  this  volunteer  in  response  to  frontal,  as  well  os  lateral 
and  oblique  impact. 

The  kinematics  of  a  linkage  selected  from  this  range  is  compared  to  the  kinematics 
observed  in  the  photography  to  evaluate  the  faith  with  which  it  reproduces  the  filmed 
data  and  the  similiarity  to  the  calculated  kinematic  trajectories  of  the  volunteer. 
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E  -  Linearized  approximat  ion  of  the  mean  square  residual  error 
Eq  -  Mean  square  residual  error  for  o  particular  geometry 

D  -  Diagonalized  matrix  of  the  second  derivatives  of  E  with  respect 

to  variations  in  the  linkage  geometry 
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X  -  Variations  in  the  linkage  geometry  according  to  the  coordinates  of  D 
-  Position  of  the  pivot  fixed  relative  to  the  first  thoracic  vertebral  body 
H  -  Position  of  the  pivot  fixed  in  the  head 
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INTRODUCTION 


The  Naval  Aerospace  Medical  Research  Laboratory,  New  Orleans,  Louisiana,  is 
currently  engaged  in  a  series  of  experiments  to  determine  human  dynamic  response  to 
impact  acceleration.  In  the  course  of  these  experiments,  human  volunteers  undergo 
short  duration  accelerations  approximating  frontal,  lateral,  and  oblique  crash  impact. 

The  resulting  motion  of  the  volunteers'  heads  and  first  thoracic  vertebral  bodies  is  monitored 
by  anatomically  mounted  clusters  of  inertial  transducers  and  photographic  targets.  It  is 
hoped  that  the  study  of  these  motions  will  contribute  to  the  design  of  safe  and  effective 
impact  protection  systems  for  use  in  all  kinds  of  vehicles. 

Perhaps  the  most  useful  result  of  such  study  will  be  the  development  of  improved 
anthropomorphic  dummies  for  use  in  testing  at  potentially  injurious  levels  of  impact. 
Differences  observed  between  the  impact  response  of  the  human  volunteers  and  that  of 
some  currently  available  anthropomorphic  dummies  show  that  the  dummy  neck  systems 
are  particularly  deficient.  (1,2)  The  subject  of  this  presentation  is  the  first  part  of  a 
two  part  approach  to  developing  a  dummy  neck  system  suitable  for  frontal,  lateral,  and 
oblique  crash  impact  testing. 

This  approach  is  developed  from  cin  earlier  attempt  to  identify  a  dummy  neck  system 
suitable  for  frontal  crash  impact  testing.  (3)  The  first  part  of  the  approach  is  a  purely 
kinematic  analysis.  All  the  data  for  each  particular  human  volunteer  is  examined  to  identi¬ 
fy  a  spatial  linkage  that  will  reproduce  all  the  head  versus  first  thoracic  vertebral  body 
positions  observed  in  the  experiments. 

There  are  three  aspects  to  this  linkage,  its  form,  its  geometry,  ond  its  articulation, 
ne  form  will  apply  to  the  whole  range  of  humanity  ond  perhaps  even  to  the  primates 
used  in  impact  testing.  The  geometry  consists  of  those  parameters  which  apply  the  general 


form  to  a  single  individual.  These  geometrical  parameters  will  vary  between  individuals 
but  will  be  constant  over  each  particular  individual's  entire  history  barring  changes  due 
to  injury,  disease,  or  advancing  age  The  articulation  consists  of  those  parameters  that 
identify  a  particular  head  versus  T-1  position.  These  articulational  parameters  vary 
continuously  with  time  and  will  become  the  variables  by  which  an  individual's  head 
versus  T-l  response  to  impact  will  be  described. 

The  second  part  of  this  approach  will  build  on  the  results  of  the  first  part.  Dynamic 
elements  will  be  incorporated  into  the  kinematic  linkages  identified  for  the  various  human 
volunteers  so  that  the  entire  assembly  will  simulate  the  dynamics  of  the  observed  impact 
response  as  well  as  the  kinematics. 

These  dynamic  elements  will  consist  of  springs,  dushpots,  and  frictions  that  will 
operate  across  the  linkage  articulations.  These  too,  will  each  have  a  general  form 

that  will  serve  across  the  range  of  experimental  subjects  but  with  precise  parameters  f 

defined  for  each  individual  subject. 

The  success  of  the  approach  depends  largely  on  the  choice  of  the  linkage  form.  A 
poorly  chosen  form  would  probably  foil  in  the  kinematic  analysis,  but  it  is  quite  possible 
that  a  candidate  form  could  satisfy  the  kinematics  and  still  make  impossible  demands  of 
the  dynamic  elements  later  on.  Since  the  choice  of  linkage  form  may  be  subject  to 

i 

reconsideration  and  this  kinematic  analysis  exercised  a  number  of  times  for  each  form, 
this  presentation  will  attempt  a  complete  treatment  of  the  kinematic  analysis  and  leave 
the  dynamic  analysis  to  the  future. 

OVERVIEW 

Since  the  kinematic  analysis  is  largely  a  study  of  the  relative  positions  of  the  head 
and  the  first  thoracic  vertebral  body;  the  primary  data  used  in  this  analysis  is  photographic. 
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The  photographic  system  used  to  acquire  this  data  is  essentially  that  described  in  ref.  4 
although  a  number  of  improvements  hove  since  been  incorporated. 

This  system  acquires  high  speed  films  of  targets  fixed  relative  to  the  volunteer's 
head  and  first  thoracic  vertebral  body  during  on  impact  experiment.  The  positions  of  the 
head  and  first  thoracic  vertebral  body  are  calculated  by  finding  those  positions  for  which 
the  theoretical  positions  of  the  target  images  in  the  film  frames  most  closely  resembles 
the  positions  actually  observed.  Given  a  kinematic  linkage  of  known  geometry,  the 
articulational  parameters  of  this  linkage  are  found  in  exactly  the  same  manner. 

Such  a  kinematic  linkage  was  selected  at  the  outset  of  this  project.  This  linkage  is 
essentially  the  two  dimensional  linkage  of  ref.  3  with  two  additional  hinges  or  pivots  as 
shown  in  figure  1  and  described  in  appendix  B. 

This  four  pivot  linkage  has,  within  the  limits  imposed  by  midsagittal  symmetry,  nme 
geometrical  parameters.  Since  the  angle  formed  by  the  axes  of  the  two  interior  pivots 
has  been  fixed  at  ninety  degrees,  only  the  eight  lengths  indicated  in  the  figure  must  be 
identified. 

The  identification  process  is  iterative.  A  likely  set  of  geometrical  parameters  is 
selected  and  compared  to  the  data  of  a  single  volunteer.  The  articulational  parameters 
are  fitted  for  each  photographically  acquired  position  of  the  volunteer  and  the  quality  of 
these  fits  is  examined  ever  the  entire  data  set.  This  examination  suggests  a  new  set  of 
geometrical  parameters  which  are  then  subjected  to  the  same  comparison.  Presumably, 
this  process  will  lead  to  the  best  set  of  geometrical  parameters  that  the  linkage  form  will 
afford  for  the  particular  volunteer. 

K  this  finol  linkage  fits  the  observed  data  well,  the  analysis  will  be  extended  to 
other  volunteers  or  even  to  the  dynamic  port  of  the  approach;  otherwise,  core  has  been 
taken  to  simplify  directing  the  analysis  to  other  forms  of  linkage. 
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The  Photographic  Data 


In  the  impact  experiments  the  volunteer  wears  two  clusters  of  photographic  targets, 
the  first  fixed  as  rigidly  as  possible  with  respect  to  the  first  thoracic  vertebral  body,  the 
second  fixed  similarly  with  respect  to  the  head  as  shown  in  figure  2.  During  impact  three 
sled  mounted  cameras  acquire  film  images  of  these  targets  at  nominal  rates  of  500  frames 
per  second. 

After  the  impact  the  film  from  the  three  cameras  is  developed  and  digitized.  In  this 
digitization  each  frame  of  film  is  examined  to  determine  the  precise  time  at  which  it  was 
exposed  and  the  position  of  each  of  the  target  images  within  the  filmplane. 

The  various  processes  by  which  the  film  images  are  related  to  the  sled  referenced 
positions  of  the  head  and  T-l  have  been  analyzed  and  quantified.  The  location  of  each 
target  in  its  cluster  is  measured  using  a  site  developed  stereovideographic  process.  The 
location  of  the  clusters  with  respect  to  the  anatomical  bodies  is  measured  similarly  using 
the  stereoradiographic  process  reported  in  ref.  5.  Optical  survey  techniques  locate  each 
camera  within  the  sled  reference  system.  Finally,  the  mechanisms  by  which  the  cameras 
produce  images  are  quantified  with  the  camera  calibration  technique  recorded  in  ref.  4. 

It  is  now  a  matter  of  arithmetic  to  set  up  algorithms  yielding  ideal  target  image 
positions  as  functions  of  the  anatomical  body  positions.  These  algorithms  may  then  be 
linearized  and  used  in  an  iterative  process  to  find  the  anatomical  positions  that  provide 
the  best  fit  of  the  ideal  image  positions  to  those  observed  in  the  digitization. 

These  anatomical  positions  can  be  solved  for  all  the  photographic  data  yielding  time 
histories  of  the  motion  of  T-l  and  the  head  in  the  sled  reference  system  in  response  to 
rash  impacts.  However,  there  are  drawbacks  to  this  direction  solution.  The  nature  of 
the  data  collection  imposes  errors  that  will  affect  this  solution  in  various  ways. 

The  most  serious  of  these  errors  are  the  small  biases  present  in  the  camera  surveys  and 
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calibrations.  These  biases  introduce  a  finite  residual  difference  between  the  observed 
image  positions  and  the  best  fit  of  the  idealized  image  positions.  Although  these  biases 
are  themselves  invariant,  the  resulting  finite  residual  can  impose  discontinuous  shifts 
in  the  calculated  trajectories  of  the  head  and  T-l . 

These  shifts  occur  because  the  motion  induced  by  the  impact  is  generally  so  large 
that  various  targets  pass  in  and  out  of  the  camera  fields  of  view  during  the  xtct  response 
The  clusters  are  arranged  so  that  sufficient  targets  are  visible  at  any  instan'  >owever  the 
set  of  visible  targets  changes  throughciut;  and,  every  time  the  target  set  ch  the  pre¬ 

sence  of  an  irreducible  residual  will  cause  a  shift  in  the  trajectory  determined  by  least 
squares  fit. 

This  phenomenon  can  be  visualized  by  invoking  the  mechanical  representation  of  a 
least  squares  fit  as  a  series  of  springs  connecting  each  observation  to  its  ideal  counterpart. 
The  best  solution  is  generally  that  for  which  the  springs  are  in  equilibrium.  If  a  spring  is 
added  to  or  taken  from  the  representation,  then  the  system  is  no  longer  in  equilibrium  and 
must  seek  some  new  equilibrium  state.  The  only  exceptions  are  those  cases  in  which  there 
i 5  perfect  agreement  between  the  observation  ord  the  ideal  counterpart  to  be  added  t©  or 
dropped  from  the  solution. 

While  these  difficulties  can  certa  nly  be  cured  by  refining  the  various  input  data  to 
eliminate  the  small  biases  that  are  the  source  of  the  problem,  the  algorithms  themseives 
can  also  be  improved.  The  solution  described  earlier  treats  the  motion  as  a  collection  of 
isolated  positions.  By  recognizing  that  the  motion  is  instead  a  sequence  of  positions,  a 
-lore  refined  algorithm  can  control  the  transitions  between  equilibrium  states  imposed  by 
ranging  target  sets. 

This  new  algorithm  shown  in  *igure  3  can  be  said  to  replace  the  springs  of  a  least 
squares  solution  with  parallel  spring  dampers.  If  on  observation  is  dropped  from  the  solution 
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the  dampers  will  impose  a  smooth  transition  to  the  new  equilibrium.  If  an  observation 
is  added  to  the  solution,  the  corresponding  spring-damper  is  set  to  the  current  equilibrium. 

This  algorithm  uses  the  straightforward  least  squares  fit  to  generate  a  series  of  small 
adjustments  for  each  of  the  photogra,:jhic  image  positions.  These  adjustments  are  weighted 
sums  of  the  residual  differences  between  ideal  and  observed  image  positions  for  prior  frames 
of  data.  These  adjustments  are  then  added  to  the  observed  image  positions  to  produce  a 
new  data  set  whose  solution  has  no  irreducible  residual. 

This  new  data  set,  christened  'at  less'  data,  is  then  itself  subjected  to  the  standard 
least  squares  algorithm.  Since  the  reiidual  differences  between  the  ideal  and  arless  image 
positions  are  almost  nonexistent,  changes  in  the  target  set  produce  no  shifts  in  the  solution. 
Furthermore,  the  solution  obtained  frxn  this  arless  data  is  asymptotic  to  that  of  the  unadjusted 
or  raw  data  approaching  it  at  a  rate  fixed  by  the  weighiing  scheme  used  to  generate  the  adjust¬ 
ments. 

The  Linkage  Analysis 

Ideally,  the  kinematic  analysis  I  an  attempt  to  infer  the  nature  of  the  constraints 
that  the  cervical  spine  imposes  on  the  relative  positions  of  the  head  and  T-l  .  To  be  sure, 
more  worldly  concerns  have  also  been  imposed,  these  being  that  the  constraints  be  mobilized 
as  some  simple  mechanism  and  further  that  this  mechanism  serve  in  conjunction  with  simple 
dynamic  elements  to  reproduce  dynam  c  as  well  as  kinematic  response;  yet,  the  kinematic 
ana'ysis  is  still  the  study  of  constraint;. 

Since  the  solution  for  head  and  T-l  position  obtained  by  the  photographic  data  system 
assumes  no  constraints  acting  between  the  head  and  T-l,  no  constrained  solution  can  obtain 
better  fit  to  the  observations.  Therefore,  the  photographic  solution  becomes  the  baseline 
for  the  study,  and  the  arless  data  obtained  by  this  solution  will  be  the  data  to  which  the 
linkages  will  be  fitted. 
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The  form  of  a  likely  linkage  candidate  is  shown  in  figure  1  and  described  in  appendix  B. 
This  is  a  four  pivot  linkage  reducing  the  head  -  T-l  system  from  12  degrees  of  freedom,  six 
each  for  the  two  anatomical  sites;  to  10  degrees  of  freedom,  six  for  say  T-l  and  four  more  to 
obtain  the  head  position. 

The  solution  for  these  10  parameters,  the  six  components  of  the  T-l  position  plus  the 
four  linkage  articulations  proceeds  in  exactly  the  same  manner  os  the  standard  least  squares 
solution  for  the  head  and  T-l  positions  except  that  instead  of  solving  for  two  independent  sets 
of  six  variables,  the  solution  obtain;  a  single  set  of  10. 

This  solution  for  the  10  parameters  will  also  yield  an  irreducible  residual  error  between 
the  constrained  idealized  image  positions  and  those  of  the  arless  data.  This  irreducible 
residual  is  due  solely  to  the  constraints  imposed  by  the  linkage  and  indicates  the  quality 
with  which  the  linkage  reproduces  the  kinematics  of  the  volunteer. 

The  purpose  of  the  linkage  analysis  is  to  find,  for  the  kinematics  of  any  single  volunteer, 
the  geometrical  parameters  with  which  the  candidate  linkage  best  approximates  the  arless  data. 
This  approximation  is  that  which  produces  the  smallest  irreducible  residual  squared  and  summed 
over  the  entire  data  set . 

This  search  for  the  geometrical  fwramefers  superimposes  a  second  least  squares  calculation 
for  the  articulational  parameters.  Bi  t,  since  orticulational  parameters  are  calculated  for 
every  position  in  the  data  set,  there  are  a  staggering  number  of  unknowns  involved.  The 
crux  of  the  matter  is  this;  The  search  for  the  geometrical  parameters  establishes  crosstc'ks 
between  the  formerly  independent  calculations  for  the  various  sets  of  orticulational  para¬ 
meters.  Instead  of  solving  n  sets  of  10  equations  in  10  unknowns,  the  task  is  now  solving  n 
times  10  plus  eight  (eight  geometrical  parameters)  equations  in  n  times  10  plus  eight  unknowns. 

Fortunately,  this  new  formulation  yields  readily  to  a  little  manipulation.  The  matrix 
arithmetic  described  in  appendix  A  obtains,  in  a  single  pass  over  the  data  set,  the  sums 
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necessary  to  solve  for  a  linearized  approximation  to  the  best  possible  set  of  the  eight  geo¬ 
metrical  parameters.  Although  the  development  indicates  that  the  articulational  parameters 
can  be  determined  by  a  second  pass  over  the  data  set,  since  the  search  for  the  geometry  is 
iterative,  this  second  pass  is  never  undertaken. 

Although  it  is  possible  to  solve  the  eight  simultaneous  equations  in  eight  unknowns 
directly  to  obtain  estimates  for  the  geometrical  parameters,  a  more  complex  approach  is 
taken.  The  main  reason  for  this  more  complex  approach  is  that  it  is  extremely  likely  that 
the  equation  set  is  very  nearly  singular.  That  is  instead  of  a  single  optimum  the  data  can 
be  satisfied  equally  well  by  a  range  o'  geometrical  parameters. 

If  such  a  range  of  solutions  is  indeed  the  case,  then  an  iteration  using  straightforward 
techniques  would  at  best  obtain  only  a  fingle  parameter  set  from  this  range  and  could  fail 
to  converge  on  an  optimum  at  all.  Thus,  the  approach  first  taken  supplied  an  interactive 
operator  with  the  results  of  the  eigenvector  analysis  described  in  appendix  C  so  that  he 
might  select  the  geometry  of  successive  iterations  directly.  After  this  technique  had  been 
exercised  a  number  of  times,  the  selection  technique  was  automated  and  the  operator  re¬ 
moved  from  the  iteration. 

The  advantages  offered  by  the  eigenvalue  analysis  is  that  it  transforms  the  eight 
simultaneous  equations  in  the  eight  geometrical  parameters  into  eight  independent  equations 
operating  each  on  a  different  iinear  combination  of  these  eight  parameters.  The  analysis 
also  shows  which  of  these  various  linear  combinations  will  most  improve  the  geometry  and 
which,  if  any,  have  little  or  no  effect  on  linkage  performance. 

The  automated  technique  solves  these  independent  equations  and  applies  the  solutions 
m  order  of  descending  importance  to  obtain  the  geometry  for  the  next  iteration.  After  the 
final  iteration,  the  investigator  obtains  a  local  optimum  linkage  and  a  statement  showing  how 
variations  of  this  linkage  will  affect  the  total  performance. 
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The  Applicotion 

This  linkage  analysis  was  applied  to  impact  data  collected  on  a  single  human  volunteer. 
This  particular  volunteer  is  a  healthy,  male,  naval  recruit  and  has  met  the  stringent  physical 
requirements  for  participation  in  the  NAMRLD  acceleration  program.  Table  1  lists  a  number 
of  anthropometric  measures  and  calculations  for  this  volunteer  as  performed  and  reported  by 
C.  Clauser  and  K.  Kennedy  of  the  6570th  Aerospace  Medical  Research  Laboratory  at  Wright 
Patterson  Air  Force  Base.  His  sitting  height,  93.5  cm,  is  in  the  67th  percentile  of  sitting 
height  for  the  population  of  naval  aviators  described  in  ref.  6. 

The  impact  data  was  taken  from  four  acceleration  experiments,  two  lateral  impacts 
(G  )  at  nominal  peak  sled  accelerations  of  6  and  11  G's  respectively;  and  two  oblique 
impacts  (G+^  x)  of  7  and  12  G's.  The  methodology  with  which  these  experiments  were 
conducted  is  essentially  that  reported  in  references  7  and  8.  In  addition  to  this  impact 
data,  three  sets  of  voluntary  motion  data  were  also  analyzed.  This  voluntary  motion  data 
is  collected  in  precisely  the  same  marner  as  the  impact  data  except  that  instead  of  under¬ 
going  impact,  the  volunteer  moves  his  head  through  a  rehearsed  trajectory.  These  trajectories 
include  pitch  motion  in  which  the  volunteer  throws  his  head  and  neck  forward  simulating  the 
purely  midsagittal  responses  observed  n  frontal  (G  )  impact;  roll  motion  in  which  the  head 
and  neck  are  rolled  from  left  to  right  while  continuing  to  face  forward;  and  yaw  motion  in 
which  the  head  is  rotated  about  the  z  axis  as  in  naying. 

The  voluntary  motion  was  incorporated  into  the  analysis  in  lieu  of  frontal  impact  data 
which  as  yet  has  not  been  collected  on  this  particular  volunteer.  The  yaw  motion  was  of 
particular  interest  because  it  was  not  observed  at  all  in  the  impact  response. 
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RESULTS 


Although  a  number  of  analyses  were  attempted  in  which  vorious  weighting  schemes 
were  applied  to  the  data  from  the  impact  and  voluntary  motion  runs,  a  single  finding  was 
common  to  all  of  them;  The  four  pivot  linkage  as  described  cannot  accommodate  all  of 
the  various  head  and  neck  motions  encountered  in  the  data.  There  is  a  conflict,  seemingly, 
between  the  voluntary  yaw  motion  and  everything  else. 

However,  since  this  yaw  motion  is  not  observed  in  any  of  the  impact  responses,  the 
linkage  is  still  considered  a  viable  candidate.  But  since  the  linkage  cannot  accommodate 
the  entire  range  of  head  and  neck  motion,  the  linkage  optimization  has  been  limited  to 
the  kinds  of  motion  observed  in  the  impact  studies. 

Accordingly,  the  yaw  motion  has  been  dropped  because  it  is  not  observed  at  all,  and 
the  voluntary  roll  motion  has  been  dropped  because  it  is  available  in  the  volunteer's  impact 
response.  The  voluntary  pitch  motion  has  been  retained  because  it  is  quite  similar  to  what 
might  be  encountered  in  Frontal  impact  and  because  there  is  as  yet  no  frontal  impact  data 
for  this  particular  volunteer. 

The  weighting  scheme  employed  took  data  evenly  distributed  throughout  each  of  the 
impact  runs  and  the  voluntary  pitch  motion  run.  The  calculation  was  such  that  oblique 
impact,  lateral  impact,  and  voluntary  pitch  were  each  weighted  equally  in  the  linkage 
optimization . 

As  expected,  no  single  optimum  was  found.  Instead,  o  range  of  linkage  suit  the  obser¬ 
vations  equally  well.  A  single  linkage  in  this  range  hos  been  mobilized  to  generate  the 
eigenvector  analysis  and  comparisons  with  the  actual  kinematics  of  the  volunteer. 

This  particular  linkage  shown  in  figure  4,  was  selected  largely  by  its  mechanical 
suitability  for  use  in  dummy  design.  However,  since  this  analysis  is  purely  kinematic  and 
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no  dynamic  work  has  yet  been  attempted,  this  particular  linkage  should  be  taken  as 
representative  of  the  performance  possible  wi  th  a  four  pivot  mechanism  rather  than  the 
last  word  in  such  mechanisms. 

The  performance  of  this  linkage  versus  each  of  the  impact  and  voluntary  motion 
experiments  is  shown  in  table  2.  The  voluntary  yaw  result  differs  dramatically  from  ‘he 
other  six  experiments.  However,  since  this  residual  is  the  root  mean  square  residual 
taken  over  the  entire  experiment,  even  this  dramatic  difference  is  understated. 

Figure  5  shows  this  rms  residual  at  various  instants  throughout  each  experiment.  But 
even  though  this  residual  is  the  criterion  by  which  the  linkage  is  evaluated,  what  these 
values  mean  in  terms  of  differences  between  the  head  and  neck  position  of  the  volunteer 
and  the  best  fit  of  the  linkage  is  not  immediately  obvious.  Figures  6  through  9  show  these 
differences  in  terms  of  magnitudes  of  translational  and  angular  displacement.  This  trans¬ 
lational  magnitude  is  given  in  meters  and  is  a  familiar  mathematical  concept.  The  angular 
magnitude,  given  in  degrees,  derives  from  the  fact  that  the  relative  orientation  between 
two  3-dimensional  coordinate  systems  may  be  expressed  as  a  single  rotation  about  a  common 
fixed  axis  as  discussed  in  references  9  and  10.  As  shown  in  these  figures,  the  yaw  motion 
leads  to  translational  and  angular  displacements  of  up  to  10  centimeters  and  30  degrees. 

The  results  of  the  eigenvector  analysis  are  summarized  in  tables  3  and  4.  This  parti¬ 
cular  treatment  attempts  to  show  how  the  linkage  performance  might  be  affected  by  varying 
certain  combinations  of  its  geometrical  parameters.  These  combinations  are  mutually  independent 
so  that  the  changes  in  linkage  performance  due  to  variations  along  a  number  of  these  combina- 
ions  is  fust  the  sum  of  the  changes  due  to  each  variation. 

Table  4  shows  just  how  large  the  range  of  acceptable  linkage  geometries  may  be.  The 
Columns  labeled  "  -»-X '  (1%)"  and  "-X1  (1%)"  show  how  the  given  linkage  can  be  modified 
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Table  2 


Run 

LX  2072 
LX  2782 
LX  2873 
LX  2904 
LX  2967 
LX  2969 
LX  2977 


Linkage  Performance  Summary 


Im 

Lateral  Impact 
Lateral  Impact 
Oblique  Impact 
Oblique  impact 
Voluntary  Pitch 
Voluntary  Roll 
Voluntary  Yaw 


Displacement 


.014  millimeters 
.0)9 
.019 
.012 
.017 
.020 
.057 


Weighted  Equo,  werghH ng  on  Mch 

average  of  thrW)  typ#f 


.017  fn  il  I  irr>efers 
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Geometrical  Parameters 
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Table  4 


Mean  Square  Residual  Versu,  Linkage  Variation, 

Elx,'L'  k.  k|2-28;  x'  1  +  E 

a  =  |  l  «  a  J  o 

(Equation  (6)  from  appendix  C) 

Eo  =  .275  x  ?<f3  MM2 


aa 

B 

X  Minima 

+  x  n°'d 

v  /  / 

(dimensionless) 

MM 

MM 

MM 

WM 

.474 

x  10"5 

.347 

x  JO"5 

0.7 

1.05 

C.3 

.677, 

.509, 

.93!  x 

« io"7 

:io-7 

IO"8 

-  .302 

-  .359 

-  .372 

x  !0"7 

x  10'7 

x  IO"6 

-  0.4 

-  0.7 

-  33.0 

6.0 

6.7 

4.1 

-c.8 

-  1.1 

_  7 

.599  x 

.863  x 

IO’8 

io"9 

. 192  x  10"6 

-  .1 53  x  !0'6 

32.0 

-  177. 

70.5 

8.8 

-  6.5 

-  36:.. 

.200  x 

kj'9 

.525  x 

io"8 

26. 

147. 

-  9* . 

.899  x 

io'10 

*  .442  x 

IO-7 

-  492. 

30. 

~  lOli. 
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Figure  4.  A  Suitable  Linkage  Candidate 
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Figure  6.  Linkage  Performance  -  T-l  Displacement 
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Figure  7.  Linkage  Performance  -  Head  Displacement 
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Figure  8.  Linkage  Performance  -  T-l  Disorientation 
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Figure  9.  Linkage  Performance  -  Head  Disorientation 


without  increasing  the  mean  square  rssidual  beyond  one  percent  of  its  present  value. 

The  first  eigenvector,  or  combination,  has  a  range  of  about  one  millimeter;  but  for 
the  eighth  eigenvector  this  range  is  well  over  o  meter. 

The  column  labeled  "X'  Minima"  is  the  theoretical  position  in  terms  of  these 
eigenvectors  of  the  absolute  optimum  linkage.  Comparing  this  column  to  other  two 
indicates  that  the  given  linkage  is  close  to  the  edge  of  the  given  range.  In  fact,  the 
figures  indicate  that  the  total  residual  could  be  reduced  by  as  much  as  25%  for  a  13 % 
reduction  in  the  rms  residual  merely  by  using  this  theoretical  solution. 

Unfortunately,  since  this  eigenvector  analysis  is  based  on  linear  approximations,  it 
should  not  be  extended  beyond  the  near  neighbors  of  the  given  linkage.  By  modifying 
the  candidate  linkage  by  five  centimeters  at  a  time  and  reinvoking  the  eigenvector 
analysis  to  reduce  the  total  residual  one  quickly  will  arrive  at  a  linkage  that  is  completely 
unworkable  and  whose  performance  is  only  marginally  better  than  that  given. 
CONCLUSIONS 

The  techniques  described  here  are  a  workable  means  of  fitting  a  proposed  linkage 
form  to  photographic  data  showing  human  kinematics  in  response  to  impact  or  in  voluntary 
motion . 

The  particular  linkage  form  selected,  a  four  pivot  mechanism,  shows  promise  for  use 
in  modeling  human  head  and  neck  response  to  frontol,  lateral,  and  oblique  impact.  Un¬ 
fortunately,  the  linkage  cannot  include  voluntary  yaw  motion  in  its  range  of  articulation. 

This  particular  inadequacy  may  well  be  ignored  as  such  yawing  motion  has  not  been 
->bs erved  in  the  impact  experiments  performed  ot  this  Laboratory. 

This  work  has  shown  that  a  range  of  linkage  geometries  will  serve  the  observed 
kinematics  equally  well .  This  finding  is  probably  generally  true  of  simple  mechanical 
simulations  of  the  human  spine  and  oervica!  column 
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The  limited  angular  motions  involved  as  well  os  the  distributed  articulation  over  a 
large  number  of  joints  should  serve  to  make  the  pivots  of  simple  models  indistinct.  This 
is  surely  the  case  in  the  two  dimensional  model  of  the  human  spine  identified  by  Mital,  et  a  l 
in  reference  11.  In  this  work  calculations  based  on  photographic  coverage  of  the  motion  of 
a  cadaveric  subject's  pelvis  and  T-l  In  response  to  frontal  impact  led  to  the  geometry  of  a 
two  pivot  linkage.  The  investigators  were  then  able  to  alter  this  geometry  considerably 
to  obtain  a  mechanically  feasible  linkage  that  still  seemed  to  serve  the  observed  kinematics. 

Since  this  work  is  the  first  part  of  a  two  part  approach,  much  of  its  utility  depends  on 
the  success  of  the  second  part,  the  dynamic  modeling  effort.  However,  it  is  absolutely 
crucial  that  this  kinematic  modeling  be  extended  to  other  human  volunteers.  Until  it 
can  be  shown  that  this  linkage  is  generally  applicable  to  the  impact  kinematics  of  human 
beings  the  work  presented  here  is  at  best  a  curiosity. 
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Ideal  value  of  the  .th  observation  in  the  nth  set  of  observations 

i  ** 


Sets  of  local  variables  that  opply  respectively  to  the  sets  of  observations 
A  set  of  global  variables  that  applies  to  oil  the  observations 


Observed  values 


Known  constants 


AAa 

ABa 

Aa 

BB 


Two  dimensional  matrices  composed  of  various  sums 


The  Kronecker  delta 


Vectors  composed  of  various  sums 
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Least  Squares  Criteria 


Find  O—i  and  <)> 
uk  m 


such  that  Z  T  (O  .  -  X  )  is  a  minimum; 

a  j  ai  a; 


that  is  find  O  , 
at< 


and  (p  that  satisfy 


i  A  A 1 1 

o 

o 

.  . 

O 

lAB,  1 '  ft' 

°lkM 

O 

|aa2I 

o 

•  • 

o 

Iab/  [2? 

02kM 

O 

o 

|aa3 

o 

,  AO  ,T  ?31 
JAB31  | 

• 

• 

» 

• 

• 

« 

* 

• 

• 

• 

• 

• 

- * —  &Q  M  1 

o 

o 

o 

1  •  • 

|4W 

.T  U  m  1 

Iabt! 

!ab2! 

|ab3I 

•  • 

|AW 

iBBl  *' 

*1  01ik  (°li  '  Xol| 


f<'*|k®2|V-V 


f '“311c  (°3(  '  Xo3|>1 


(2) 


*  » 
^  J aaMfk  (°a  Mf  -  Xoa 


v  r  I  b_.  fOn.  -  X  .)] 
r-  *•  1  Qjm  Q|  on  i 


)  ! 


I  a 


oai 


|AAal 


°ail 


'"“a1  ‘  ?  V  bai» 


|bb|  -  II  b  b 

i  a  a;rr'  ^i1 


Ibis  rnanix  formulation  is  potentially  quite  large  and  also  contains  many  zero  elements.  Solving 

or  0  and  <f)  by  straightforward  matrix  inversion  techniques  would  be  at  best  clumsy  and  m- 
Ok  m 

ficient.  At  worst,  limitations  on  computer  core  and  word  length  would  render  such  a  solution 
meoningless,if  not  impossible. 
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Although  it  is  possible  to  obtain  a  simpler  formulation  by  direct  manipulation  of  equation  (2) 
a  clearer  development  might  be  obtained  as  follows: 


Instead  of  solving  directly  for  the  best  fits  ofc^  and<J>  ,solve  for  the  best  fit  of 
0  ^  as  o  function 
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where  the  superscript  *  indicates  a  best  fit 

Now  substitute  this  expression  forO^  (<|>^)  back  into  equation  (I)  and  rearranae  to  aet 
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Now  find  the  best  fit  for  <J> 


-1 


(5) 


(Z((6  -A  )  (6  -A  )b  (0_  -  x  )  j ) 

1  '  vu  Ovu  v fj  Ovw  uur  Uw  OQw 


The  unique  structure  of  the  matrices  A  makes  them  symmetric  and  idempotent,  that  is 
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As  a  result  6-  A^  is  also  idempotent  so  that  equation  (5)  may  be  rewritten 
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Note:  In  general,  repetition  of  indices  implies  summation  except  that  index  represented  by  a 
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X  Vector  in  three  dimensional  coordinates 

<D  Transformation  matrix  or  pivot  rotation 

C  Transformation  matrix 

Superscript  primes  denote  successive  coordinate  systems 
Superscript  T  denotes  matrix  transpose 
Subscript  o  denotes  pivot  location 
Subscript  1,2,3  identifies  pivot 
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The  Linkage  Candidate 

The  constraint  that  a  pivot  put;,  on  the  relative  positions  of  two  coordinate  systems 
may  be  written 

X  -  0  (X  -  X  )  +  X  ( 
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where  0  can  be  described  os 
=  a 


Cos  0  0  Sin  0 
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C  =  C  0  c 


-Sin  0  0  Cos  0 


That  is  rotation  about  a  single  axis  oriented  by  C  in  the  X^  system  and  by  C^  in  the 
X  system.  X  and  X  represent  ct  single  point  on  this  axis  in  X^and  X_  respectively. 

A  chain  of  coordinate  systems  constrained  by  a  sequence  of  pivots  may  be  described 
for  two  pivots  as 

*  *,2-*o.>  *  *o>  (3) 

and  for  many  pivots  as 

X  =  0  ,  (0  ,(...©  <Xn  -  X"  )+Xn_1  ....)  +  x‘  ,-x'  )  +  X  (4) 
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Since  the  interior  coordinate  systems  of  these  linkages  are  of  little  interest,  they 
may  be  oriented  and  located  as  necessary  to  simplify  these  expressions.  Rewriting  equation 
(3)  yields 

£>,.£,  £2^2 h  <-x"--V 
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=.1  =yl  =1  -a2  -al  -al 
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C^  Cj  is  a  coordinate  rotation  and  may  be  described  as  successive  rotations  about  y, 
carried  'z',  and  carried 'y'  yielding 


(6) 
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I  I 

furthermore,  the  vector  (X  -  X  ,)  moy  be  similorly  broken  down  to  components  porallel 

-  o2  -  ot 

to'y'j,  'z',  ond^.  Since  the'yj  component  is  unaffected  by  rotations  about  y^  it  may  be 

passed  backward  through  C  and  0  ..  The'yl  component  may  be  passed  forward  through 

=  yl  =  y  I  l 

I 

C  _  and  0  *  in  a  similar  manner. 
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I  ^^1 

If  we  now  select  the  X  system  so  that  z  is  perpendicular  to  the  two  pivot  axes 

I  » 

contained  in  this  system  then  C  ^  and  may  be  reduced  to  identify  matrixes  and  the 

)  I  l 

remaining  component  of  (X  -  -  X  )  lies  along  the  z  axis.  Equation  (5)  becomes 
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vectors  in  the  X  and  X  systems  respectively  they  may  re  incorporated  into  X  and  X 

—  a  1  -  a  2 

respectively  to  yield 
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If  bilateral  symmetry  and  pivo  axes  parallel  to/y','y4' ,  and  each  other  are  imposed 
then  the  still  general  expression  (8)  oecomes 


This  statement  has  five  parameters:  X  Z  Z  ,  X  ,  and  Z  which  correspond 

at  at  or  or 

T]  ,  L,  and  in  figure  1. 

The  four  pivot  expression  resui "ing  from  equation  (4)  is 
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The  interior  systems  and  X^'s  tray  be  redefined  to  obtain 
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This  derivation  began  between  the  two  interior  pivot  axes  and  proceeded  oit  ward 
to  the  ends  of  the  linkage.  Since  components  were  passed  outward  only,  two  of  the  interior 

I 

vectors  analogous  to  0,  0,  Z  in  equation  (7)  must  also  have  xaxnponents. 

This  expression  may  be  simplified  further  by  imposing  the  following  requirements: 
3i'ateral  symmetry,  the  first  axis  is  parallel  fo'y',  the  second  is  perpendicular  fo'y,  the 
third  axis  is  perpendicular  to  y  and  the  fourth  is  parallel  to  y  .It  is  also  assumed 
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that  for  the  case  of  zero  rotation  about  the  second  and  third  axes  that  y  is  parallel  to 

l  l  1 1 

'y'  .  The  expression  now  becomes  i 
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The  quantity  Z  is  an  interesting  casualty.  The  requirements  on  the  axis  orientations 

and  bilateral  symmetry  combine  to  force  it  to  zero. 

The  rotation  C,  satisfies  the  perpendicularity  requirement.  It  may  be  dropped  from 

=J. 

the  expression  by  altering  the  forms  of  the  interior  rotations  to  rotations  about 'x'and'z' axes 
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in  this  work  C  is  further  assumed  to  be  an  identity  matrix  reducing  the  description  of  this 

-y 

I  I  III  til  I  It  l  till 

system  to  eight  fixed  parameters.  These  are  X^,  Z^,  X  ,  Z  ,  X  ,2  ,  X^,  and  Z^ 

which  correspond  to  T^,  T  ,  L^,  L H^,  and  in  figure  1. 


B-5 


SYMBOL  TABLE  -  APPENDIX  C 

A  Symmetric  matrix 

B  Vector 

£  Seal  or  -  Mean  Square  Error 

*  Vector  -  Variation  in  linkage  geometry 

Orthogonal  transformation 
^  Diagonalized  matrix 
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APPENDIX  C 


Eigenvalue  Analysis 

The  essence  of  the  linkage  analysis  is  that  it  seeks  to  obtain  that  linkage  geometry 
for  which  a  single  sum  is  a  minimum.  This  sum  is  the  total  of  the  squared  residua  I  differ¬ 
ences  between  the  arless  data  and  the  best  fit  to  this  data  by  the  articulational  parameters. 

As  can  be  seen,  there  are  two  best  fit  procedures  going  on  simultaneously;  one 
involving  the  geometrical  parameters  and  the  other  the  articulational  parameters.  These 
become  the  global  and  local  variables  respectively  of  appendix  A. 

The  analysis  depends  on  obtaining  linearized  approximations  of  the  functional 
relationships  of  these  residual  differences  with  the  various  parameters.  These  linearized 
relationships  can  then  be  used  to  infer  the  behavior  of  that  single  sum  of  squared  residuals 
for  a  range  of  geometries. 

Of  course  the  inference  here  depends  on  the  reliability  of  the  linearized  approxi¬ 
mations  of  the  residuals  versus  the  parameters,  and  these  are  reliable  only  in  the  vicinity 
of  the  parameters  for  which  the  approximations  were  derived. 

Fortunately,  the  situation  is  not  quite  circular,  the  approximations  may  be  derived 
for  some  parameter  set  chosen  largely  by  intuition.  Then  the  linea  inference  can  be  trusted 
to  identify  a  better  parameter  set  that  is  still  in  the  vicinity  of  the  initial  guess.  Successive 
derations  of  this  process  may  then  lead  to  on  optimum  linkage  for  which  the  sum  of  the 
squared  residuals  is  truly  a  minimum. 

The  inference  is  in  the  form  of  a  single  matrix  equation  called  a  quadratic  form. 

A,.  X,  X.  -  2  B,X.  +  E  =  E  (X)  0) 

lj  I  j  1  i  o 

vVhere 

x  represents  the  difference  between  any  geometry  and  the  set  of  parameters  currently  considered: 
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A  Is  an  eight  by  eight  symmetrical  matrix  and  B  on  eight  vector  all  of  whose  components 
are  calculated  in  the  various  sums  of  appendix  A;  is  the  squared  and  summed  residuals 
for  the  particular  parameters  being  considered;  and  E(X)  is  the  proiected  value  for  thot 
sum  at  a  geometry  removed  by  X  from  the  current  parameters. 

The  value  of  X  for  which  E(X)  is  a  minimum  can  be  derived  by  applying  a  little 
differential  calculus 


i 


2  A..  X.  -  2B. 
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(2) 


A..X.  B. 
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The  direct  solution  would  then  be: 

X.-A.."1  B. 
i  <1  I 


(3) 

(4) 


which  presumably  would  be  added  to  the  current  set  of  geometry  to  repeat  the  analysis 
around  this  projected  minima. 

If  the  projected  value  of  X  is  large,  it  is  quite  possible  that  the  actual  value  of  E 

will  be  larger  than  the  projected  value  and  in  fact  even  larger  than  E  itself.  It  is  very 

o 

likely  in  such  a  case,  that  repeated  calculation  would  diverge  from  the  best  geometry. 
Furthermore,  if  the  equations  are  nearly  singular,  that  is  if  the  matrix  A  is  ill  conditioned, 
X  will  almost  certainly  be  large. 

These  difficulties  can  be  circumvented  by  applying  an  eigenvector  analysis.  This 
technique  conjures  up  on  eight  dimensional  countour  mop  in  which  each  of  the  eight  di- 
n  ensions  corresponds  to  one  of  the  geometrical  parameters  and  the  contours  correspond  to 
the  values  of  E(X).  These  contours  are  concentric  ellipsoids  whose  corresponding  values 


of  E(X)  diminish  toward  the  center. 

The  eigenvector  analysis  yields  the  shape  of  these  ellipsoids  and  their  orientation 
in  space  by  diagonalizing  the  matrix  A. 

A..  =  C, .  D  C.. 

ij  ki  kl  I  j 

where  D.  ,  =  0  k  / 1  (5) 

kl 

and  C  represents  an  affine  transformation . 

The  matrix  C  represents  a  transformation  to  a  new  set  of  coordinates  in  the  eight  dimensional 
space.  The  axes  of  this  new  system  are  parallel  to  the  axes  of  the  concentric  ellipsoids. 

The  off-diagonal  elements  of  D  are  equal  to  zero,  but  the  diagonal  elements  are  inversely 
proportional  to  square  of  the  corresponding  semi  axes  of  the  ellipsoids. 


Applying  this  transformation  to  equation  (I)  yields 
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Where  the  prime  indicates  that  the  vectors  hove  been  rotated  to  the  new  coordinates.  This 
new  equation  reduces  the  contribution  of  changes  in  linkage  geometry  to  the  total  residual 
error  to  eight  independent  contributions. 

The  advantage  of  this  new  formulation  is  that  these  eight  independent  contributions 
may  not  be  equally  important.  A  small  variation  along  one  of  these  primed  axes  may  produce 
a  large  change  in  E  while  even  larae  variations  along  other  axes  may  have  only  the  most 
negligible  effect. 

I  ♦ 

The  X  for  which  EfX  )  is  a  minimum  is 

I  I 

X  -  B  /D  no  sum  on  a  17) 

a  a  a  a 

The  corresponding  decrease  in  E  for  each  component  of  X  is 

’  2 

-  E  =  fB  )  /D  no  sum  on  (81 

a  a  aa 
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and  the  second  derivative  of  E(X  )  is 

&  =  2Dfla  0  -  a 

dx'adx'0  0  0  /o 


(9) 


The  eigenvector  analysis  reduces  the  eight  dimensional  search  to  eight  independent 
one  dimensional  searches.  The  results  of  these  eight  searches  may  be  Invoked  selectively 
at  the  whim  of  the  investigator  within  the  range  of  the  linear  approximation. 

Furthermore,  since  Dq  is  directly  proportional  to  the  second  derivatives  of  the 
eight  solutions,  it  serves  to  rank  them  in  order  of  importance.  The  solution  is  best  defined 
along  the  axes  for  which  DQa  is  greatest.  As  Dao  approaches  zero,  displacement  along  the 
corresponding  axis  becomes  less  significant  indicating  that  a  range  of  geometries  may  serve 
equally  well.  Although  negative  values  of  Daa  are  theoretically  impossible,  round  off 
errors  and  the  like  may  result  in  infinitesimal  Daa‘s  being  calculated  as  slightly  less  than 
zero.  Such  Daa  may  be  Ignored. 

At  the  outset  of  the  linkage  study  the  results  of  the  eigenvector  analysis  were  displayed 
to  the  investigator  so  that  he  might  select  the  next  candidate  geometry  directly.  On  -he 
basis  of  this  experience  the  routines  were  restructured  to  invoke  the  eight  eigenvector  solu¬ 
tions  in  order  of  descending  Daa  until  a  maximum  displacement  was  reached  or  a  maximum 
decrease  in  expected  E  was  exceeded.  This  scheme  worked  well  in  practice. 
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